On the nuclear dependence of the /i — ► e 
conversion branching ratio 
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Abstract 



The variation of the coherent branching ratio R^ e (ratio of the /i - — > e~ reaction 
rate divided by the total muon-capture rate) through the periodic table is studied 
by using exact muon wave functions. It was found that, by using very heavy nuclei 
(e.g. 197 Au, the SINDRUM II target) as [i~ — > e~ conversion stopping-targets, the 
above ratio is favored by a factor of about four to five than by using light ones (e.g. 
48 Ti, chosen as PRIME target). 
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The present status of the exotic neutrinoless /i — > e conversion in nuclei, 



fi + {A,Z)^e- + (A,zy, (1) 



has comprehensively been discussed recently [1-3] both from experimental [4- 
9] and theoretical [10-15] point of view. Process (1) violates the Li, (i = [/,, e) 
quantum numbers and has been proposed [12] as one of the best probes to test 
the existence of charged-lepton flavor conservation, among a great number of 
similar processes predicted by modern gauge and supersymmetric theories. 

Several experiments have been designed to explore process (1) and performed 
at PSI (and earlier at TRIUMF) on 48 Ti, 208 Pb and 197 Au targets [1,4,7,8]. 
They have put so far, only bounds on the branching ratio R^. Presently, the 
published best upper limit is R T ^ < 6.1 x 1CT 13 [7]. The ongoing SINDRUM 
II experiment (PSI) is now using 197 Au as stopping target and the extracted 
preliminary limit constitutes an improvement over the previous one set on a 
heavy target ( 208 Pb [4]) by two orders of magnitude [7,8]. The planned MECO 
experiment (Brookhaven) is going to use 27 A1 target and a very intense pulsing 
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muon-beam [5,6] to reach a sensitivity of roughly R™ < 2 x 1CT 17 [5,6] which 
implies an improvement of the present limits by about three orders of mag- 
nitude. It should be mentioned that very recently a new fj,~ — > e~ conversion 
experiment on 48 Ti (PRIME) was announced to be performed at KEK [2,9] 
aiming to push the limit down to R^ l e < 1CT 18 . In all these experiments, the 
signature of reaction (1) is a single electron with energy E e = m M — e& — m e 
(neglecting recoil) where e& is the muon binding energy in the Is orbit and 
(m e ) the muon (electron) mass. 

From the theoretical point of view [12-17], process (1) constitutes a very 
good interplay between atomic, nuclear, particle and non-standard physics 
[14,15]. The fj,~ — > e~ conversion Hamiltonians which result in the context 
of many extensions of the standard model proposed up to now, in general, 
give rise to coherent and incoherent processes [16]. In several models, like 
those for which the isoscalar couplings of the vector and scalar interactions 
are not very small [3,16], the coherent action of all nucleons in \x~ — > e~ 
leads to enhancement of conversion electrons thereby making it potentially a 
sensible indicator for lepton flavor violation (LFV) effects [10-13]. In addition, 
the g.s. — > g.s. transitions are favored due to Pauli blocking effects which 
prevent the formation of excited states. The present work is motivated from 
the necessity to investigate the nuclear physics aspects of process (1). We 
study the nuclear structure dependence of the coherent branching ratio R^ e 
throughout the periodic table by performing exact calculations of the muon- 
nucleus overlap integrals. 

The expression for R^ e , to leading order in the non-relativistic reduction, for 
the coherent process has been written in the form [16] 

_ F, e (A,Z) G% Pe E e \M^ 



where T Me stands for the fjT — > e~ conversion rate and T MC for the total rate 
of the ordinary muon-capture, ^ + (A, Z) — > + (A, Z — 1) [18,19]. The 
factor G\j1 corresponds to non-photonic mechanisms and for photonic ones 
it should be replaced by the ratio (47ra) 2 /g 4 , where a is the fine structure 
constant. In Eq. (2), p e denotes the outgoing-electron momentum connected 
to the excitation energy E x of the daughter nucleus (E x = Ej — E gs ) through 
the relation 

p e w q = m ^- tb - E x . (3) 



q = |q| is the magnitude of the momentum transfer (we neglect the electron 
mass). The quantity Q of Eq. (2) depends very weakly on the nuclear structure 
and, in principle, it contains scalar (S), vector (V), axial- vector (A), pseudo- 
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scalar (P), and tensor (T) coupling terms [16]. Especially for photonic dia- 
grams, which are the main concern of this work, Q is rather nuclear-structure 
independent. Hence, the main nuclear physics aspects of R^ e are accumulated 
in the last fraction of Eq. (2). 

The matrix elements Ai^\ a = V, S, A, P,T, which enter the expression of 
R^ e are defined by 

= (flEQlitie-^HvjM, (4) 

r = for isoscalar and r = 1 for isovector operators, respectively, where \i) 
the initial and |/) the final nuclear state, ^(ij) represents the muon wave 
function evaluated at the position of the j th target-nucleon. The functions 
contain the spin-isospin dependence of the p~ — > e~ operator [3]. For 
the coherent process (|/) = in the case of scalar and vector interactions, 
My \ are written in terms of the ground-state proton, neutron densities as 

M^ s (q) = J [ Pp (r) ± p n (r)] e'^^r = F p (q) ± F n {q) , (5) 

the (+) sign corresponds to r = and the (-) to r = 1 channel, where 

F P M = J PpM e " iq r ®»(r)d 3 r ■ (6) 

The proton (neutron) density p p (p n ) is normalized to the atomic number 
Z (neutron number N) of the nucleus in question. For our purposes here, the 
required densities p p are taken from experiment [20]. For photonic mechanisms 
only protons of the target-nucleus contribute and hence, M.y\ = M.y\ = 
F p (q). 

For light and medium nuclei (see discussion of the results below) can 
be reliably evaluated in a straightforward way by factorizing outside the in- 
tegrals of Eq. (4) a suitably averaged muon wave function ($^). Under these 
conditions Eq. (4) is approximated by 

m { : ] = (<£>;*></ 1 Ee;o>-^ | <> = (&;)m£), (7) 

where involve the pure nuclear physics aspects of the pT — > e~ pro- 

cess. Equation (7) for photonic diagrams, is written as My'ig = F P {q) = 
ZF p (q). For the mean muon wave function (<&* s ), a simplified expression 
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(see e.g. Ref. [18] and Eqs. (22), (23) of Ref. [17]) was used in muon capture 
studies by many authors [18]. In previous estimations of the branching ratio 
R^ e , the same expression for ( < &^ s ) was adopted [13] for both the numerator 
and the denominator of Eq. (2) in order to reduce the uncertainties inserted 
via the use of Eq. (7) in /i~ — > e~ and /i~ — > processes. 

In the special case of g.s. — > g.s. transitions for spin zero (J=0) light nuclei, 
Myg are determined by the elastic scattering (monopole) nuclear form factors 

F p , F n [13] and they are simply given byMg(g) = ZF p (q)±NF n (q). In 
the general case of nuclei with ground-state spin J ^ (e.g. 27 Al, the MECO 
target, with J = |), My^ s contain, in addition to monopole (L = 0) form 
factors, contributions arising from other multipoles (L = 2,4,...). However, 
the latter contributions are not significant [3]. 

The main goal of the present work, was to study systematically the exact 
nuclear-structure dependence of R^ e (A, Z) by computing the integrals of Eq. 
(5) for the dominant coherent — > e~ conversion and investigate the in- 
fluence on R^ e of the following effects: (i) the approximate evaluation of the 
muon-nucleus overlap integrals [see Eq. (7)] and (ii) the neglect of the muon- 
binding energy e h in Eq. (3). The latter assumption implies that for all nuclei 
in the coherent mode we have 

q = p e ^ m^/c = 0.534 fm' 1 , E e « 105.6 MeV . (8) 

In this work, the rather simple photonic mechanism, where only the target- 
protons contribute to R^e, was examined. The main steps followed in the 
calculational procedure and the results obtained are briefly discussed below. 

In the first step, the variation of the quantity Z\Fz(q)\ 2 through the periodic 
table was studied. As is well known, Weinberg and Feinberg [10] using the 
above approximations have shown that, for the coherent /i~ — > e~ rate it 
holds 

oc Z\F NN {q)\ 2 . (9) 

In Fig. 1(a) the results obtained for Eq. (9) in the following cases are illus- 
trated: 

(i) In the first case, experimental form factors F z (q) = F NN (q) [20] were used 
at the values of q given by Eq. (8) (dashed-dotted line) and by Eq. (3) (dotted 
line). 

(ii) In the second case, the phenomenological expression for F NN (q) [see Eq. 
(16) of Ref. [10]] was used in Eq. (9) (solid line). 
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The common feature of the three curves in Fig. 1(a) is the fact that they show 
a maximum in the region of A pa 60 (copper region) as had been estimated 
in Ref. [10]. This behavior was generally adopted by the experimentalists ex- 
ploring the /i _ — > e~ process [1,2]. As it can be seen, in the region of light 
nuclei, the three curves nearly coincide, which means that the approximation 
of Eq. (8) is reasonable in this region, but for medium and heavy nuclei where 
6b becomes significant, the obtained rates when q is given by Eq. (3) are larger 
by about a factor of two than those when q is given by Eq. (8). 

In the second step of the calculations, the variation of the quantity R^ e / Q p h 
[see Eq. (2)], which contains the main nuclear-structure dependence of R^ e , 
was studied. Two cases were distinguished: 

(i) In the first case, the mean muon wave function was inserted in Eq. (2). 
This gives 

^.Mrf^ ^W (10) 



The results obtained in this way, are presented in Fig. 1(b) and correspond 
to the two choices of coherent momentum transfer discussed before: (a) by 
using the values of Eq. (8) (solid line) and (b) by using the values of q given 
by Eq. (3) (E x = for coherent mode) (dashed line). In choice (a) we see 
that R^ e (A, Z) presents a maximum at the region of A pa 130 which is not in 
accordance with the estimation of Ref. [10]. This is due to the fact that in Ref. 
[10] the gross nuclear dependence of r^ c was considered as linear in Z which 
is equivalent to a constant Primakoff function fcp{A, Z). It is well known, 
however, that fop(A,Z) is strongly dependent on the mass excess (see e.g. 
Ref. [13]). In this work, in order to minimize such uncertainties, experimental 
data for the total muon capture rates T MC were used [19]. From the results 
of choice (b) we see that R^A, Z) shows a slow increase up to the heaviest 
nuclei. 

(ii) In the second case, the explicit muon-nucleus overlap integral of Eq. (4) 
was used (numerical integration) in Eq. (2). Then, the variation of R^(A, Z) 
relies on the expression 

^eOcW^^ (11) 
Q J- nr. 



This requires the use of the exact muon wave function calculated as in Ref. 
[21]. Here the values of q are provided by Eq. (3). The results are represented 
by stars (*) in Fig. 1(b). We see that, the exact evaluation of the muon- 
nucleus overlap integrals of Eq. (5), shows linear increase of R^ e as function 
of A (or Z [21]). For heavy and very heavy nuclei (region of 197 Au and 208 Pb), 
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the use of the exact muon wave function in Eq. (11), gives much larger rates 
than the approximation of the averaged muon wave function Eq. (9). The 
comparison is much worse if e& is neglected in the kinematics. One must notice 
that, the presence of the factor g~ 4 in Eq. (11) (photonic mechanism), inserts 
an additional (A,Z) dependence on the branching ratio R^ e which had been 
previously overlooked. 

In conclusion, the gross variation of R^ e obtained by the exact results, shows 
a linear rise with A which can be attributed to the coherent effect. Due to this 
behaviour of R^ e (A, Z), very heavy nuclei are favored by a factor of about five 
to be used as — > e~ conversion targets. It is worth mentioning that, the 
above factor may, in some cases, be partly compensated by other experimental 
advantages of some specific muon-stopping targets [5]. We also remark that in 
the present work we neglected relativistic atomic effects related to fi~ — > e~ 
process [14] which might influence a bit (depending on the nucleus) the results 
for R^. 

In summary, we have investigated the dependence of the fj,~ — > conversion 
branching ratios on the nuclear parameters A and Z throughout the periodic 
table. This exotic process is an interesting and important one to be stud- 
ied, since stringent bounds for the charged-lepton flavor violating parameters 
already exist and significant improvements over these limits are feasible in 
the not-too-distant future from the SINDRUM II, MECO and PRIME experi- 
ments. We performed direct calculations of the muon-nucleus overlap integrals 
(for photonic mechanisms) and found that R^ e {A,Z) keeps increasing up to 
the very heavy nuclei. This shows that, by using such isotopes (e.g. 197 Au, 
present SINDRUM II target) as targets in fj,~ — > e~ conversion experiments, 
R^ e is favored by a factor of about four to five than by using light isotopes. 
Finally, the present exact results were exploited to test the validity of previous 
approximations made on the study of the /i — e conversion and shed more light 
in this situation. 

The author wishes to thank Dr. Y. Kuno for financial support and Dr. Andries 
van der Schaaf for fruitful discussions on SINDRUM II experiments. 
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FIGURE CAPTION 



Fig. 1. Variation of the fi~ — > e~ conversion branching ratio, i? Me , through 
the periodic table assuming that the nuclear-structure dependence of R^ e is 
described: (a) by Eq. (9) [10] and (b) by Eqs. (10) and (11). For details see 
the text. 
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